Azimuthal flow of decay photons in relativistic nuclear collisions 



O 
O 

(N 



m 

(N 
> 

o 
in 
o 
\o 
o 



o 

=5 



X 



Biswanath Layek, Rupa Chatterjee, and Dinesh K. Srivastava 
Variable Energy Cyclotron Centre, 1/AF Bidhan Nagar, Kolkata 700 064, India 

(Dated: February 9, 2008) 

An overwhelming fraction of photons from relativistic heavy ion collisions has its origin in the 
decay of 7r° and n mesons. We calculate the azimuthal asymmetry of the decay photons for several 

azimuthally asymmetric pion distributions. We find that the kr dependence of the elliptic flow 

o 

parameter V2 for the decay photons closely follows the elliptic flow parameter v 2 evaluated at 
p T rs kr + 8, where S tn 0.1 - 0.2 GeV, for typical pion distributions measured in nucleus-nucleus 
collisions at relativistic energies. Similar results are obtained for photons from the 2-7 decay of n 
mesons. Assuming that the flow of n° is similar to those for 7v + and tv for which independent 
measurements would be generally available, this ansatz can help in identifying additional sources 
for photons. Taken along with quark number scaling suggested by the recombination model, it 
may help to estimate V2 of the parton distributions in terms of azimuthal asymmetry of the decay 
photons at large kr- 



I. INTRODUCTION 

The azimuthal flow of particles produced in relativistic 
heavy ion collisions has provided a strong evidence for 
the creation of a hot and dense system very early in non- 
central collisions . The importance of this observation 
stems from the fact that pressure gradients generated in 
the system at very early times, transform the eccentricity 
in co-ordinate space for such collisions to the momentum 
space for distribution of the produced particles . The 
Pt dependence of the elliptic flow parameter vi at lower 
transverse momenta has been quantitatively explained 
using hydrodynamics 0, The observed saturation of V2 
at higher pt has been attributed to effects of viscosity jj| 
or incomplete thermalization 5]. The observed scaling 
of vi with the number of valence quarks in the hadrons 
is understood in terms of the recombination model for 
hadronization 0,0- 

It has recently been suggested that the study of the 
elliptic flow of thermal photons @ may provide valuable 
insight into the build-up of azimuthal flow with time. 
This will require a subtraction of the contribution of pho- 
tons from the decay of tt° and r\ mesons produced in the 
collisions 

The deviation of elliptic flow of inclusive photons from 
that of decay photons, will thus confirm the presence of 
additional sources of photons. We shall see later that 
the elliptic flow of photons from the decay of tt° and r\ 
mesons may also provide useful estimates for V2 of the 
partons in the frame-work of the recombination model. 

In the present work, we suggest an ansatz for the eval- 
uation of the elliptic flow of photons from the decay of 
pions, which could be useful for such studies. In the 
next section we study the transverse momentum depen- 
dence of the elliptic flow parameter for decay photons 
for several momentum distribution functions for the ir° 
and r\ mesons. In section III; we evaluate the t>2 for pho- 
tons arising from 2-7 decay of 7r° and r\ mesons, which 
are formed from recombination of partons to get a di- 
rect relation with V2 for partons. Finally, we give our 
conclusions. 



II. PHOTONS FROM DECAY OF tt° AND n 
MESONS 

Consider a ir° having four-momentum p and mass m 
decaying into two photons. The momentum distribution 
of photons in an invariant form is given by [To| : 



(1) 



where k is the four-momentum of the photon. 

Thus, the Lorentz invariant cross-section for photon 
production is 
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Let us write the pion distribution as: 



S(p-k- \m 2 ) . (2) 
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for a non-central collision of identical nuclei, where 
i> 2 (pt) is the momentum dependent elliptic flow pa- 
rameter. We consider two typical parameterizations for 
the pion distribution: an exponential distribution and a 
power-law distribution. For the first case, we assume 



da 



PT dpT dy 



exp(-yJp 2 T + m 2 /T a ) x cxp (-y 2 /2a) (4) 



where the slope parameter Tq is of the order of 290 
MeV and the width parameter of the rapidity distri- 
bution a is taken as ~ 4. For the power-law distribution, 
we take [HI 
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(5) 



where po ~ 5 GeV and n is about 29. The elliptic flow 
parameter V2 is approximated as 



v 2 {pt) = a [1 - exp(-p t /6)] 



(6) 
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FIG. 1: (Color on-line) Upper Panel: Spectrum of photons 
from the decay of 7r° for an exponential (Eq. and a power- 
law (EqU^J distribution for pions. Lower Panel: Elliptic flow of 
photons from decay of n° having elliptic flow given by Eq. © 
at y = 0. f(pr) stands for the momentum distribution of 
the 7T°. The symbols give the result of numerical calculation, 
while the curves give the fits. 



which approximately reflects the increase of v 2 with px 
for lower transverse momenta and saturation of its value 
at larger pt- We have taken a = 0.2 and b — 1 GeV for 
the first set of calculations. All the results given in the 
present work are for photon rapidity equal to zero. 

Eq.|21is then evaluated numerically and ^(^t) for the 
resulting photon distribution obtained. In Fig. ^ we give 
our results for the spectrum of photons from decay of 
7T° along with the variation of transverse momentum de- 
pendence of their elliptic flow parameter. We find that 
for both distributions, the resulting «2(&t) for the decay 
photons can be approximated as: 
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FIG. 2: (Color on-line) Same as Fig. Qfor r\ mesons. 



and S w 0.1 - 0.2 GeV, for kr > 0.2 GeV, to an accuracy 
of better than 1 -3 %. 

This result can be understood as follows. The S- 
function in Eq. ^ provides that 

m T k T cosh(y p ~y k )-p T k T cos(<t) p -(f)k)-^m 2 = 0, (9) 

in an obvious notation. Next, we note that once the 
momentum of the pions is large, the opening angle for the 
decay-photons will be very small [To| . Thus, for example, 
in the extreme case, we have, y p — ~ and <j) p — 
4>k ~ or 7r, for the photons which are co-linear with the 
pion. This leads to the solutions: » px and kx ~ 
(m 2 /4pT)(l — m 2 /4p T ) for large pr- 

However, in general, the two photons between them 
will cover the entire range of the allowed transverse mo- 
mentum. Both the photons, which will be almost co- 
linear with the pion, will "inherit" the v 2 of the pion, 
which has a larger transverse momentum. The photon 
with the lower kx will, however, be submerged in a much 
larger yield of photons coming from decay of tt° having 
lower px , as the transverse momentum distribution of pi- 
ons falls steeply with increase in px- Thus, in general the 
v 2 of a photon with a given transverse momentum will 
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be larger than the vi of a pion with the same transverse 
momentum. This accounts for the shift, Eq. |SJ 

This argument will hold even if the tj 2 (px) for the pions 
decreases with increase in px, which is likely for larger 
px- The only difference would be that tj 2 for the decay 
photons will now be smaller than the V2 for the pions at 
the same transverse momentum (see later; Fig.[5J|. 

The spectra and the transverse momentum dependence 
of the elliptic flow parameter for the 77 mesons are given in 
Fig. [3 The rich structure seen for the rj mesons at smaller 
transverse momenta has its origin in the large mass of the 
77 mesons. Thus, we find that the x>i for decay photons 
coming from the 2-7 decay of 77 mesons, approaches the 
V2 for the mesons only at larger transverse momenta (> 
0.8 - 1.0 GeV). For the simple parameterizations of the 
momenta considered here, t> 2 is negative for very low kx 
as the decay photons are distributed away from the ma- 
jor axis in the momentum space when the px of the 77- 
mesons is low and the opening angles are large. 

Of-course, a more complete analysis would include the 
photons coming from the 3-7r° decay of the 77 mesons, 
which will mainly populate the low kx window. While 
a simulation of the contribution of this process to the V2 
can be straight-forward using standard event generators, 
a direct evaluation will not be easy as it would involve a 
10-dimensional numerical integration. 

However, we can easily make an estimate of the range 
of the transverse momenta of photons, where the photons 
arising from decay 77 — > 3ir° — > 67 would contribute. The 
maximum kinetic energy of a 7r° in the rest frame of the 
r] meson is given by |l4j: 
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4t7C 
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80McV. 



(10) 



Thus, for example, the maximum energy of a pion 
emerging from an 77 meson having a momentum of 2 GeV 
would be of the order of 0.62 GeV and would contribute 
to photons having momenta less than that. Next, we re- 
call that 77/77° is of the order of 0.44 for px > 2 GeV and 
is expected to drop to zero at px = 0, at least accord- 
ing to calculations based on PYTHI A 15]. This coupled 
with the branching ratio of about 0.3 for this mode of de- 
cay should ensure that the contribution of photons from 
the 3-7T° decay of 77 mesons would be limited to smaller 
momenta. A more detailed study of this effect would 
definitely be useful. 

Let us return to the discussion of the V2 of photons 
from the 27 - decay of ir° and 77 mesons. We now plot the 
ratio of the elliptic flow parameters of the calculated de- 
cay photons and the mesons as a function of their trans- 
verse momenta in Fig. [3] We see once again that V2 
for the two are quite similar when transverse momenta 
exceeds 1 GeV. Similar results for the transverse momen- 
tum dependence of v 2 /v% had earlier been noted by the 
WA98 group ^(|. We now understand this in terms of 
the momentum shift <5 (Eq. [SJl • 

Considering the likely usefulness of the relation given 
by Eq. [H] we have examined it to determine the depen- 
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FIG. 3: Transverse momentum dependence of the ratio of 
elliptic flow parameters for decay photons and the n° or 77 
meson. 
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FIG. 4: Variation of the shift parameter S with the slope pa- 
rameter (To) and the power-law parameter n for the photons 
coming from the decay of ty°. 



dence of 5 on the slope-parameter To or the power-law 
parameter n (see Fig0J. We see that, in general, the shift 
5 is smaller if the distribution is steeper. This may be 
considered as a generalization of the Stcrnheimcr 's 
prescription for azimuthally asymmetric distribution of 
pions and their decay photons. 

We now take a more general shape of the transverse 
momentum dependence of the tj 2 for pions, which rises 
with px, reaches a maximum, and then starts decreasing. 
We illustrate our discussion (see FigEJ) by taking, 



V2(pt) = 0Ap T exp(-/j T ) 



(11) 



We again see that the transverse momentum dependence 
of V2 for the decay photons is obtained from the V2 for 
the pions by a shift of about 0.27 GeV. 

It has been recently reported that the so-called 
quark-number scaling of the elliptic flow parameter sug- 
gested by the recombination model works much more ac- 
curately when plotted in terms of "transverse kinetic en- 
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FIG. 5: Elliptic flow parameter for photons from the 2-7 de- 
cay of 7T°. The symbols denote the calculated vi while the 
dashed curve denotes the fit. The «2(pr) for pions is taken 
to rise first, reach a maximum and drop for larger transverse 
momenta. 
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FIG. 7: (Color on-line) Elliptic flow parameter for photons 
(solid curve) from the decay of 7r° (dashed curve) obtained 
using the recombination model. 



collisions is simply related (Eq|HJ to the elliptic flow pa- 
rameter for 7T°. As we can get independent information 
about the later in terms of the flow of ir + and ir~ , these 
results could be useful in identifying additional sources 
of photons @- 



III. RECOMBINATION MODEL AND DECAY 
PHOTONS 

As remarked earlier, one of the more interesting re- 
sults of elliptic flow studies has been the approximate 
quark- number scaling of the hadronic V2- This finds a 
natural explanation in terms of the recombination model 
for hadronization of the quark gluon plasma 0, Q , which 
suggests that in the region of pr where the recombina- 
tion of the partons dominates the process of hadroniza- 
tion, that is pr < 4 (6) GeV for mesons (baryons) and 
the effect of mass is minimal, Vi obeys a simple scaling 
law: 



ergy" or tut — m, where m is the mass of the hadron. In 
order to see the consequences of this behaviour of the el- 
liptic flow parameter, we now replace pr with nix — m in 
Eq.|SJ Next, taking the exponential distribution function 
(Eq. for transverse momenta, we calculate the elliptic 
flow parameter for decay photons from tt° and r\ mesons 
(see Fig.©. 

Identifying fc-r of the photons with the so-called "trans- 
verse kinetic energy" , we again see that V2 for photons is 
simply related to v-i for the pions by a small shift of w 0.1 
GeV in the transverse kinetic energy of 7r° . We also find 
that we can easily describe the «2 of the photons from 
the decay of r\ mesons by a small negative shift. 

Thus, we see that the elliptic flow parameter of decay 
photons from 7r° which make an overwhelming contri- 
bution to the inclusive photons in relativistic heavy ion 



v 2 (pt) ~ riv\{p T I 'n) (12) 

where n is the number of valence quarks in the hadron 
and v\ is the elliptic flow parameter for them. 

The shifted-scaling of the i>2 of the decay photons with 
the t>2 for the ir° and rj seen earlier for a variety of mo- 
mentum distributions and momentum-dependences of V2, 
taken along with the valence-quark scaling of the v 2 for 
hadrons, holds out the promise of the "lowly" decay pho- 
tons providing a direct measure of the v\ of the partons. 
In order to examine this possibility, we now explicitly 
calculate the v 2 for photons from the decay of ir° and ry, 
which in turn, are formed by the recombination of par- 
tons. 

The transverse momentum distribution of mesons 
formed by the recombination from thermal parton dis- 
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FIG. 8: (Color on-line) Elliptic flow parameter for photons 
(solid curve) from the 2-7 decay of ^-mesons (dashed curve) 
obtained using the recombination model. 



tribution can be written as Q: 
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In the above mr is the transverse mass of the hadron, 
m q =260 MeV is the quark mass, T = 175 MeV is the 
transition temperature, and r]T is the transverse rapidity 
such that tanh?/T = 0.55, appropriate for a system cre- 
ated in collision of gold nuclei at ^/sjvn =200 GeV. While 
writing the above we have utilized the (5-function for the 
asymptotic form of the perturbative meson distribution 
amplitude. 

We now return to the present calculation. In these 
illustrative studies, we account for the elliptic flow of 
the mesons by multiplying the above distribution with 
[1 + 4i>2(pt/2) cos(20)]. We have calculated v^ipr), the 
elliptic flow for the partons, using Eq. 81 of Fries et al. 0. 

In Fig.[7| we give our results for the elliptic flow param- 
eter for pions obtained, using the recombination model 
along with the V2 for the photons coming from the decay. 
We should not take the results for the vi for pions for 
lower px too literally, as the pion mass is much smaller 
than 2m q which goes into making it. However, for larger 
Pt this difference is not very relevant and we can trust 



the model. We see that v\ closely follows the vl{ as be- 
fore. 

The mass of rj mesons on the other hands is large 
enough to be free from this complication. We see that 
once again, the decay photons closely follow the az- 
imuthal asymmetry of the rj mesons at larger pt (see 
Fig. 0. 

IV. SUMMARY 

We have studied the azimuthal asymmetry of distribu- 
tion of photons originating from the 2-7 decay of ir° and rj 
mesons produced in relativistic heavy ion collisions. Sev- 
eral distribution functions and parameterizations of the 
transverse momentum dependence of the elliptic flow pa- 
rameter, including those coming from the recombination 
model have been considered. 

We have empirically found that the elliptic flow pa- 
rameter for the photons from the decay of pions is quite 
close to the same for the pions evaluated at a transverse 
momentum shifted by about 0.1 - 0.2 GeV. This has its 
origin in the small opening angle for the photons for pi- 
ons having large momenta. Similar results are found for 
photons from the 2-7 decay of rj mesons, at pt > 1 GeV. 

These empirical findings could be useful in identifying 
additional sources of photons, as the elliptical flow of 
7T should be similar to that for tt + and tt~ . The flow 
of rj mesons could perhaps be approximated with those 
for kaons based on considerations of their masses and 
the number of valence quarks. In any case, in general, 
the yield of eta mesons is much less compared to that 
for pions. Thus the presence of an additional source of 
photons will be indicated by the deviation of the vi for 
inclusive photons from the V2 of decay photons from the 
2-7 decay of ir° and rj mesons, which could be obtained 
using the V2 for charged pions and kaons. 

We add that a relation 

El 

between the pt integrated 
anisotropy for ir° and photons from its decay, obtained 
using simulations, has been successfully used to analyze 
the azimuthal anisotropy of decay photons in the WA98 
experiment |19|. 

Thus, an experimental verification of the empirical 
findings for the transverse momentum dependence of the 
elliptic flow reported in the present work, may be of con- 
siderable interest. 

In general, such a study would also account for the 
acceptance and efficiency of the photon detectors, which 
is beyond the scope of this work. 
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